Zero permeability and zero permittivity band gaps in ID 
metamaterial photonic crystals 

Ricardo A. Depine jj] and Maria L. Martinez- Ricci * 
Grupo de Electromagnetismo Aplicado, Departamento de Fisica, 
Facultad de Ciencias Exactas y Naturales, 
Universidad de Buenos Aires, Ciudad Universitaria, 
Pabellon I, C1428EHA Buenos Aires, Argentina 

Juan A. Monsoriu § 

Departamento de Fisica Aplicada, Universidad Politecnica de Valencia, 4-6022 Valencia, Spain 

Enrique Silvestre " and Pedro Andres ** 
Departamento de Optica, Universidad de Valencia, 46100 Burjassot, Spain 

(Dated: February 2, 2008) 

Abstract 

We consider layered heterostructures combining ordinary positive index materials and dispersive 
metamaterials. We show that these structures can exhibit a new type of photonic gap near frequencies 
where either the magnetic permeability u or the electric permittivity e of the metamaterial change 
signs. Although the interface of a medium with zero refractive index (a condition attained either when 
[i = or when e = 0) is known to give full reflectivity for all incident polarizations, here we show 
that a gap corresponding to u = occurs only for TE polarized waves, whereas a gap corresponding 
to e = occurs only for TM polarized waves. These band gaps are scale-length invariant and very 
robust against disorder, although they may disappear for the particular case of propagation along the 
stratification direction. 
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Photonic band gap (PBG) materials allow electromagnetic field propagation in certain fre- 
quency bands, but not in others, their essential feature being the periodic arrangement of 
high-contrast electromagnetic properties. The simplest version of a PBG material is the ID 
planar layered stack, which has been commonly used in optics in the form of filters or Bragg 
reflectors 

Almost every available natural material has been used to construct PBG structures. Possibil- 
ities have been widened with the recent advent of metamaterials (MMs) , artificially constructed 
composites exhibiting electromagnetic properties that are difficult or impossible to achieve with 
conventional, naturally occurring materials . Key representatives of this new class of materi- 
als are MMs with negative index of refraction, a property that arises in media with a negative 
electricpermittivity together with a negative magnetic permeability in the same frequency 
range Multilayered PBG structures containing MMs have attracted much attention in re- 
cent years . Two new kinds of PBG, fundamentally different from the usual Bragg 
gaps originating from interference in the periodic structure, have been identified: the "zero 

n n 

averaged refractive index" gap [5] and the " zero effective phase" gap [6( . The first one appears 
in multilayers combining ordinary materials (positive refractive index) and MMs with negative 
refractive index, whereas the second one appears in multilayers containing two different single- 
negative (permittivity- or permeability-negative) MMs. In contrast to Bragg PBGs, both the 
"zero averaged refractive index" and the "zero effective phase" gaps remain invariant to scale- 
length changes and are robust against disorder. PBGs with these characteristics can be useful 
as photonic barriers in quantum well applications, such as multiple channeled filtering 

Since Bragg PBGs result from interference between waves reflected at different layers, mate- 
rial dispersion does not play a crucial role in their appearence. In other words, Bragg PBGs in 
ID multilayers always exist, even when the dispersion of the constituent materials is neglected. 
On the other hand, the existence of a "zero averaged refractive index" PBG highly relies on 
the dispersive characteristics of the material layers, since for nondispersive materials this PBG 
would cover all frequencies, except for singular points [5]. In this Letter we show that MM mul- 
tilayers can also exhibit a new kind of PBG, where material dispersion -and not interference- 
plays the key role. These new PBGs correspond to frequency regions where the constitutive 
parameters of the MM change their signs. Since they are not based on the Bragg interference 
mechanism, they are also invariant to scaling and even more insensitive to disorder than the 
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"zero averaged refractive index" and the "zero effective phase" gaps. 



If losses are neglected, the refractive index of a MM n(uS) = J fi(oj)e(oj) is zero at the values of 
the frequency u where the magnetic permeability fi(uj) or the electric permittivity e{uj) change 
their sign. Therefore, at these frequency values one could expect to find gaps for both TE and 
TM polarization, since it is well known that the boundary of a zero refractive index material 
exhibits full reflectivity 0, 0]. In contrast to this expectation, we show that at frequencies 
where \l — a gap occurs for TE polarized waves, but not necesarily for TM polarized waves, 
whereas at frequencies where e = a gap occurs for TM polarized waves, but not necessarily 
for TE polarized waves. Moreover, we show that these gaps can be completely absent for the 

particular case of propagation along the stratification direction, which could explain why its 

rn 

presence went unnoticed in previous studies (see, for example, Fig. 2 in [5|). 

We consider the ID periodic structure created by layers of two different materials: a conven- 
tional dielectric with permeability permittivity e\ and thicknesses d\, and a MM with per- 
meability /i2(<^), permittivity 62(00) and thicknesses d,2- The stratification direction is the y axis 
and we consider wave propagation in the x — y plane. Let the function f(x, y) represent the z- 
directed component of the electric field for the TE-polarization case (electric field parallel to the 
layers), and the z-directed component of the magnetic field for the TM-polarization case (mag- 
netic field parallel to the layers). The propagating waves in the periodic structure have the form 
of Bloch modes, for which the fields satisfy the condition f(x, y + d) = /(0, y) expi(k x x + Kd), 
where d = d\ + c?2 is the period of the structure, k x is the wave vector component along the 
layers and K is the Bloch wave number. For two-layered periodic structures, the dispersion 
relation K(u, k x ) can be found explicitly from 

cos(Kd) = £ = cos(kiydi) cos( y k2yd 2 ) — 
1 



®2k\y <J\k%y 



sin^iydi) sm(k2yd 2 ) , (1) 



where the index j = 1,2 indicates the layer, <jj = fij for TE polarization or <jj = ej for TM 
polarization, kj y = kj — k%, and kj = uorij/c are wave numbers in each media with refractive 
indexes %. The quantity £, half the trace of the matrix characterizing the unit cell translation 
operator determines the band structure. It takes real values for lossless media and real k x . 
Regimes where |£| < 1 correspond to real K and thus to propagating Bloch waves. In regimes 
where |£| > 1, K has an imaginary part, therefore the Bloch wave is evanescent, and this regime 
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corresponds to forbidden bands (or gaps) of the periodic medium. The band edges are those 
regimes where |£| = 1. 

Let us consider frequency regions where the refractive index n 2 (oj) goes through zero, that 
is, where a constitutive parameter of the MM -either /i 2 (cu) or e 2 (o;)- changes its sign. Note 
that the term between square brackets in the right hand side of eq. (JIJ becomes singular under 
these conditions. For propagation along the stratification direction (k x = 0), k 2y — > and 
the singularity is compensated. To be specific, and taking into account that \i 2 and e 2 do not 
generally change their sign at the same frequency value, we assume that e 2 — » but fi 2 7^ 0. In 
this limit, k 2y — > and the quantity £ for both TE and TM polarizations adopts the form 

£ « cos(fcidi) - -— a — // 2 ^2 sin(^ 1 di) , (2) 
2 c y \i\ 

which can take any real value, either outside or inside the interval [—1, 1]. A similar conclusion 
can be obtained when fi 2 —>■ but e 2 7^ 0. Thus, gaps near regions where a constitutive 
parameter changes its sign do not occur automatically at normal incidence. For example, the 
multilayer considered to obtain Fig. 2 in Ref. Q], exhibits Bragg gaps and a "zero averaged 
refractive index" gap, but there is no gap in the regions where fi 2 and e 2 change signs. This 



can be observed in the band diagrams shown in 
following geometric and constitutive parameters 



igs. lb and If, obtained for k x = and the 



12 mm, g? 2 = 6 mm, €\ — /ii — 1, 



<*(/) = 1 + 9 2 _ f 2 + 1L5 2 _ f 2 , 
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where / is the frequency measured in GHz. 

For oblique incidences, on the other hand, k x ^ 0, k 2y — > ±ik x 7^ and the second term 
in the right hand side of eq. 1 remains singular when /i 2 = (TE polarization) or when 
e 2 = (TM polarization), thus originating new gaps, as shown in Fig. 1 for 9 = 30° and 60° 
(k x = k\ sin#). The non singular behavior of £ at normal incidence and its singular behavior 
for oblique incidences is illustrated in Fig. |21 for both polarizations. We conclude that these 
new gaps always appear for non zero propagation angles 9 in layer 1: the /i 2 = gap appears 
for TE, but not for TM, polarization; whereas the e 2 = gap appears for TM, but not for TE, 
polarization. 



The projected band structure corresponding to this example is shown in Fig. |H1 White 
regions indicate forbidden bands where there are no electromagnetic modes, regardless of K. 
The width of the "zero averaged refractive index" gap does not change appreciably for TE 
polarization, but it gets narrower when 9 — > 90° for TM polarization. On the other hand, the 
fi 2 = and the e 2 = gaps increase with the angle 9. 

Like the "zero averaged refractive index" gap p and the "zero effective phase" gap [6], the 
gap arising from /i 2 = or e 2 = differs fundamentally from the usual Bragg gap. First, 
its central frequency is independent of the lattice constant, while all Bragg gap frequencies 
must scale with the lattice constant. This fact is shown in Fig. HJ where we compare the 
transmittance through 16 unit cells corresponding to the periodic multilayer considered in Fig. 
El for 9 = 45°, with the transmittance obtained for similar structures with lattice constants 
scaled by factors of 2/3 and 4/3. While the Bragg gap shifts upward or downward in frequency, 
the "zero averaged refractive index" gap -near 2.3 GHz (TE) and 3.55 GHz (TM)- and the 
new gaps -near 3.55 GHz (TE) and 3.9 GHz (TM)- remain unchanged. Second, the Bragg 
gap, relying on interference mechanisms, is destroyed by strong disorder, while the /i 2 = and 
the e 2 = gaps, relying on constitutive properties of the MM, should be expected to be robust 
against disorder. In Fig. EJwe compare the transmittance through 16 unit cells corresponding 
to the periodic multilayer considered in Fig. 01 for 9 = 45° and to similar structures with 
thickness fluctuation of ±3 mm and ±6 mm, each ensemble averaged over 24 realizations. As 
expected, the Bragg gap is destroyed by disorder, but the "zero permeability" and the "zero 
permittivity" gaps survive. 

To summarize, we have shown that periodic multilayers containing MMs can exhibit a 
new type of photonic gap near frequencies where either the magnetic permeability /i or the 
electric permittivity e of the MM change signs. In contrast to Bragg PBGs, originating from 
interference mechanisms, it is the material dispersion of the MM what mainly determines 
the appearence of these new gaps. Therefore, like the "zero averaged refractive index" 
and the "zero effective phase" gaps, the new gaps also remain invariant to scaling and 
insensitive to disorder. Although electromagnetic waves would be completely reflected at the 
boundary of a medium with zero refractive index, our results show that multilayers with MM 
constituents do not automatically exhibit PBGs at frequencies where the refractive index is 
zero: for propagation along the stratification direction, zero refractive index PBGs can be 
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completely absent, whereas for oblique propagation they may emerge for different polarizations, 
depending on which of the constitutive parameters /i(a>) or e(u) makes the refractive index zero. 
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FIG. 1: Band structure for TE and TM polarizations and different angles of incidence corresponding 
to a periodic stack with air layers = ei = 1, d\ = 12 mm) and MM layers (//2 and £2 given by 
eq. El d>2 = 6 mm). The left column shows the frequency behavior of the constitutive parameters 
corresponding to layer 2. 




FIG. 2: Half of the trace of the unit cell translation matrix, £, as given by eq. ^ for the structure 
considered in Fig. ^ 9 = 0°, 30°, and 60°. The left column corresponds to TE polarization, whereas 
the right column corresponds to TM polarization. 
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FIG. 3: Projected band structure for the MM multilayer considered in Fig. ^ 




FIG. 4: (color online) TE and TM transmittance through 16 unit cells corresponding to the band 
structure in Fig. [31 for 9 = 45° (solid line), the same structure but the lattice constant scaled by 2/3 
(dashed line), and the same structure but the lattice constant scaled by 4/3 (dotted line). 
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FIG. 5: (color online) TE and TM transmittance through 16 unit cells corresponding to the perfectly 
periodic multilayer considered in Fig. [31 for 9 = 45° (solid line), to a similar structure with thickness 
fluctuation of ±3 mm (dashed line), and thickness fluctuation of ±6 mm (dotted line), averaged over 
24 samples. 
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